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ABSTRACT 

It is believed that planets are formed by aggregation of dust particles suspended in 
the turbulent gas forming accretion disks around developing stars. We describe a mech- 
anism, termed 'Stokes trapping', by which turbulence limits the growth of aggregates 
of dust particles, so that their Stokes number (defined as the ratio of the damping time 
of the particles to the Kolmogorov dissipation timescale) remains close to unity. We 
discuss possible mechanisms for avoiding this barrier to further growth. None of these 
is found to be satisfactory and we introduce a new theory which does not involve the 
growth of small clusters of dust grains. 

Subject headings: Planet formation, accretion, turbulence, caustics. 



1. Background 

It is widely believed that planets are formed by aggregation of dust particles in an accretion 
disk surrounding a growing star; the fact that solar planets have orbits which are roughly circular 
and coplanar with the Sun's equator is r eadily explained by this inodel. According to this picture, 
planet formation is a two-stage process (ISafranovlll969l : lArmitagd 120071 ) . The final stage must be 
driven by gravitational forces, but initially the density of dust particles is not sufficiently high for 
gravitational instability to overcome centrifugal forces. It is believed that the dust particles initially 
aggregate by random collisions, and that when the aggregates have grown to a sufficient size, they 
sink to the mid-plane of the accretion disk. Gravitational collapse to form planets can commence if 
the density at the mid-plane of the disk becomes sufficiently high. In this paper we are concerned 
with the first ('kinetic') phase of planet formation in this standard model. 



The kinet i c agg r egation process i n protoplanetary disks (reviewed by iDominik et al\ (|2007l ) ; 
Henning et al\ (|2006l ) : lArmitagd (|200?1 )) has distinctive features. The dust particles are surmised to 
be in a highly turbulent environment (at least during periods when the star is accreting material at a 
substantial rate), because laminar fiows could not dissipate energy sufficiently rapidly to account for 
observed accretion rates. Unless some mechanism enables the particles to fuse together chemically, 
they must form very weak aggregates (bound by van der Waals and possibly also electrostatic forces) 
until gravitational collapse is well advanced: this implies that the aggregates are very fragile unless 
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they can be heated to high temperatures. Also, the aggregates of dust particles may be fractal 
structures, with a fractal dimension D. 

The principal purpose of this paper is to argue that the standard model for planet formation by 
aggregation of dust grains is highly problematic, because the dust grains are unable to aggregate in 
the turbulent environment of the accretion disc. Furthermore, we argue that the turbulence would 
prevent even very large aggregates from settling to the mid-plane of the disc. The difficulties are 
so severe that it is necessary to develop an alternative theory. 

A note about terminology is in place here. We describe the condensed material in the inter- 
stellar medium as consisting of dust grains, which may form aggregates, bound by van der Waals or 
electrostatic forces. We argue that these aggregates may c/wsier together in space, without coming 
into contact. We also describe properties of small fragments of solid material in a turbulent gas in 
a more general context, without needing to specify whether they are dust grains or aggregates: in 
such cases we refer to particles suspended in the gas. 



2. Summary of results 



This paper uses recently acquired insights ( Wilkinson et a/.ll2006 : Falkovich et a/.|[2002 ; Wilkinson k, Mehlig 



20051 ) into the dynamics of turbulent aerosols to draw conclusions about the aggregation of dust 
particles in a standard model for a protostellar accretion system. An important parameter is the 
Stokes number, St, which is the ratio of the time scale for an aggregate of dust grains to be slowed 
by drag forces (the 'stopping time') to the dissipative correlation time of the turbulent forcing. 
Note that our definitio n of the Stokes number differ s froin that used in so me other works on planet 
formation (for example. lYoudinl (j2003l ): lBauer et al\ (j2007l ): lYoudinl (|2007l )). which define the Stokes 
number as the ratio of the stopping time to the orbital peri od, although ou r defin i tion is the one 



which i s used in t h e flui d dynamics literature, for example in lFalkovich et al\ (j2002l ) : iDuncan et al. 
(I2OO5I ): bee et al\ (I2OO6I 1. 



Our principal conclusions may be summarised as follows. We argue that the relative velocity 
of the dust aggregates, and therefore their rate of collision, suddenly increases by several orders of 
magnitude when aggregates grow to a size such that their Stokes number is of order unity. One 
consequence of this effect is the formation of a 'Stokes trap': aggregates with Stokes numbers 
significantly larger than unity will collide with sufficiently high relative velocity that they fragment 
upon collision, thereby replenishing the population of aggregates with smaller Stokes number. 
This process would result is a stable population of aggregates with a Stokes number which has 
a distribution spanning St ~ 1. The production of planets therefore depends upon the particle 
aggregation process being able to escape from this Stokes trap. We consider four scenarios by 
which escape from the Stokes trap might happen: 



1. Clustering. Particles in a turbulent flow with Stokes numbers of order unity have been 
shown to cluster together. This clustering effect has been ascribed to particles being 'centrifuged' 
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away from vortices (lMaxevlll987ll. More recent wo r k has shown that th e particles cluster onto a 
fractal set (jSommerer k, Ottlll993l : iBec et a/.ll2006l : iDuncan et al\ l2005l ) of dimension Dd- Note 
that there are two fractal dimensions in this problem: the dust particle aggregates may be fractals 
of dimension D, and these aggregates cluster onto a set of dimension Dc\. It is conceivable that this 
clustering effect might create a means of escaping the Stokes trap, but we have not been able to 
identify a mechanism which would be sufficiently effective in the model which we have considered. 

2. Temperature variations. The temperature of the dust aggregates could vary dramatically 
as they are swept back and forth between the surface and the mid plane of the accretion disk by 
turbulent motion of the gas. High temperatures can facilitate chemical reactions, as can the con- 
densation of water and other volatile materials when the dust aggregate encounters low temperature 
regions. Chemical processes could cause the dust aggregates to become sufficiently tightly bound 
that that are less easily fragmented in collisions. Also, condensation tends to occur more readily 
in corners, favouring the formation of compact aggregates. Below we argue that the accretion disk 
is sufficiently optically thick at the crucial stage of its development that there can be a marked 
difference between the surface and mid-plane temperatures. However, unless the star is accreting 
at a high rate, at distances significantly greater than 1 AU the temperatures are low everywhere 
across the profile of the disk. We conclude that chemical processes are probably not significant 
unless accretion is occurring at a very high rate. 

3. Fractal aggregates. The Stokes trap can also be avoided if the dust particles form aggregates 
with a fractal dimension D less than two: in this case we find that the Stokes number doe s not 
increase as the size of the aggregate increases. Some experiments on dust grains (jBlum et a/.ll2000. : 



Krause &; BlumI 12004 ) have shown evidence that they form aggregates with a fractal dimension 
approximately equal to 1.4, others (jWurm &: BlumlllQDg ) obtained values close to 1.9. It should be 
noted that these experiments only observe the dust particle aggregates for a short time after their 
creation, whereas in protoplanetary processes the aggregates have plenty of time to relax to a more 
compact state, which would be energetically favoured due to increased contact area. The structure 
of dust aggregates in real protoplanetary systems is therefore highly uncertain. The hypothesis that 
the dust aggregates are tenuous structures with D < 2 does not appear to provide a satisfactory 
resolution, because these tenuous aggregates will always be advected with the gas, and will not 
sink to the midplane. We also show that tenuous aggregates may be torn apart by shearing forces. 

4. Quiescence. It might be expected that the relative velocities of the dust particle aggregates 
are reduced if the rate of accretion is reduced. We write the rate of accretion onto the star as 
M = 10~^AMyr~^, where M is the solar mass, and the coefficient is chosen so that A = 1 
corresponds to a typical rate of accretion (lArmitagd 120071 ). We investigate the scaling of the 
relative velocity of suspended particles with the scaling parameter A, and find that it increases as 
A decreases. We conclude that the Stokes trap is still present in relatively quiescent disks, where 
A is small. 



Our calculations are based upon estimates using Kolmogorov's theory for turbulence (|Frisch 
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19971 ) in combination with a standard model for an accretion disk (jShakura &: Sunvaevlll973l ). In 
our application of the Kolmogorov theory we assume that the size scale of the largest eddies, L is 
proportional to the sc ale height of the disk, H, writing L = iH (we show that this is equivalent 
to the 'a-prescription' (jShakura &: Sunvaevlll973l ')). We also consider how the results depend upon 
the additional scaling parameter i. 

Because of the large uncertainties in the parameters of the model, our calculations are to be 
interpreted as illustrations of the Stokes trapping principle, rather than as quantitative predictions. 
Order of magnitude changes in the uncertain parameters will not alter our principle conclusion 
about the effectiveness of the Stokes trap in limiting dust aggregation. For simplicity we confine 
attention to the case where the heating of the disk is active (that is dominated by dissipation in 
the disk rather than by direct illumination from the star). Our calculations lead to the estimate 
that typical particle sizes of Stokes-trapped aggregates at 10 AU from the star are of the order of 
25/xm. 

In principle it is possible to gain information about the sizes of particles suspended in accre- 
tion disks around yo ung stars from spectroscopic studies, but this is an ill-conditioned problem 
([Carciofi et a/.ll2004l ). with uncertain spectroscopic data. Studies of the 'silicate feature' (at ap- 
proximately lO^m wavelen gth) in a survey of spectra of accretion disks around young stars by 
Kessler-Silacci et al\ (120071 ) are interpreted as indicating typical particle sizes bet ween 3 and lOAtm 



betwe en 10 ^ AU and lOAU from the star. Similar estimates were obtained by Ivan Boekel et al 



(j2003l ) and lApai et all (j2005l ). By contrast, spectral analysis in the mm-range (probing predomi- 
nantly the outer regio ns of the disk) has been interpreted as indicating that particles larger than 



50/xm may be present (jWood et al. 



2002 



Bauer et a/.ll2007l ^. 



Even if the particle growth process can circumvent the Stokes trap, the particle aggregates must 
collapse to the mid-plane of the disk for gravitational instability to become effective. If we make the 
assumption that the integral scale of the turbulence is comparable to the thickness of the accretion 
disk (that is, if £ is of order one), we show that advection of particles by the largest turbulent eddies 
presents another very significant barrier to planet formation: turbulence will prevent aggregates 
of centimetre size from settling to the mid-plane, unless the rate of accretion is very substantially 
smaller than 10~'''Myr~^. 

We conclude that there is no satisfactory theory for the formation of planets by aggregation of 
sub-micron sized dust grains. An alternative theory is therefore required. In the concluding section 
we introduce an alternative hypothesis. 
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3. Collision processes in turbulent aerosols 



Throughout we assume that the drag force on a particle is proportional to the difference in 
velocity between the particle and the surrounding gas, so that the equation of motion is 



'y{u{r,t) 



(1) 



where r is the position of the particle and u{r, t) is the fluid velocity field (until the particles come 
into contact). This equation is familiar in the context of Stokes's law for the drag on a sphere, 
where the damping rate 7 is proportional to the kinematic viscosity 1^. 

We assume that t he gas has a multi-scale turbulent flow, which is described by the Kolmogorov 
theory of turbulence (jFrischl 119971 ) . The length scale and time scale associated with the smallest 
eddies where energy is dissipated are denoted by rj and r respectively. According to the Kolmogorov 
theory of turbulence, these quantities are determined by the rate of dissipation per unit mass, £, 
and the kinematic viscosity v. 



1/4 



1/2 



(2) 



If the turbulent motion is driven by forces acting on a length scale L ^ rj, the v elocity fluctu ations 
of the fluid have a power-law spectrum for wavenumbers between 1/L and 1/r/ (|Frischlll997l l. 



We define the Stokes number to be St = l/(7r). Very small particles are advected with the 
gas flow, but in addition they have small random velocity due to Brownian diffusion, which can 
cause them to colli de at a rate TZd ■ They ca n also collide due to the shearing motion of the flow, an 
effect discussed by Saffman fc Turner (jl956l ) in the context of rainfall from turbulent clouds. These 
advective collisions occur at a rate TZ^ which turns out to be negligible compared to Brownian 
diffusion in the present context (see tables [3] and H]) . 

When the Stokes number approaches unity, there is a dramatic increase in the relative velocity 
of suspended particles due to the formation of fold caustics in their velocity field, illustrated for the 
case of one spatial dimension in figure [H When faster particles overtake slower ones, the manifold 
representing the phase-space distribution of the particles develops folds, and the velocity field of the 
particles goes from being single valued to multi- valued (three- valued, in this illustration). Before 
the caustics form the relative velocity of the particles is due to their Brownian diffusion. After the 
folds have formed the relative velocity may be orders of magnitude higher. The rate of collision TZ 



for a s uspended particle in a turbulent flow is approximated by a formula presented bv lMehlig et al 

n = nd + n^ + exp(-^/st)7^g . (3) 



Here TZg is the collision rate predicted by a model introduced bv I AbrahamsonI (|l975l ). often termed 
the 'gas- kinetic model', in which the suspended particles move with velocities which become uncor re- 



lated with each other and with the gas flow. A less precise formula was suggested bv lFalkovich et al 



(|2002l ). and the idea that there is a dramatic change in the relative velocity of particles at St ~ 1 
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Fig. 1. — Illustration of fold caustics. Shown are results from numerical simulations in one spatial 
dimension of the equation of motion Xi = 'y{u{xi,t) — Xi), where u{x,t) is a random velocity field 
which is a smooth function of x and t. Here Xi denotes the position of particle number i. The top 
panel shows particle trajectories Xi{t) as a function of time t. The bottom panels show a different 
representation of the dynamics: the particle velocities x = v{x,t) versus their positions x at three 
different times. The velocity v{x, t) of the suspended particles is initially single- valued, but becomes 
triple-valued in the region between two fold caustics, which are singularities of the projection of 
the phase-space manifold onto coordinate space. In the region between the caustics the suspended 
particles have a relative velocity and their trajectories cross. 



appears to have been originally proposed bv iMarkiewicz et a\ (|l99lh . The exponential term in ([3]) 
describes the fraction of the coordinate space for which the velocity field is multi-valued: A is a 
'universal' dimensi onless constant. The non-ari alytical dependence of the rate of caustic production 
on St was noted by IWilkinson Mehlid (120051 ^ and recent simulations of Navier-Stokes turbulence 
suggest that A a 2 ( Falkovich Sz Pumiil 20od ) . The rate TZg greatly exceeds TZd, but the gas-kinetic 
theory is only applicable when the velocity field of the suspended particles is multi-valued. The 
exponential term arises because the formation of caustics is determined by a process involving es- 
cape from an attractor by a diffusion process, similar to the Kramers model for a chemical reaction. 
The exponential term is therefore analogous to t he Arrhenius term expi—E /kf^T) in the expression 
for the rate of an activated chemical reaction (jWilkinson &: Mehlig |2005| ) . The abrupt increase 
of the collision rate as the Stokes number exceeds a threshold value was first noted in numerical 
experiments bv ISundaram &: Collins! (|l997l ). 
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At large Stokes number the rate of collision of particles of radius a is 

7^g = A-Ko^niAv) 



(4) 



where n is the number density of particles and Av is the relative speed of two suspended particles 
at the same position in space (and angular brackets denote averages throughout). In a multi scale 
turbulent flow, when 7r <C 1 the motion of the suspended particles is underdamped relative to the 
motion on the dissipative scale, but it is overdamped relative to slower long-wavelength motions 
in the fluid. The relative velocity of two nearby particles is a result of the different histories of 
the particles. If we follow the particles far back in time to when they had a large separation, 
their velocities were ve ry different, b u t thes e velocity differences are damped out when the particles 
approach each other. iMehlig et al\ (120071 ) surmised the variance of the relative velocities using 
the Kolmogorov scaling principle. In the following we consider, for simplicity, only the case of a 
symmetric collision, where both particles have the same damping rate, 7 (in the case of collisions 
between particles of very unequal sizes, we may still use the estimate below, using the smaller of 
the two values of 7). 

When the particles are underdamped relative to the smallest dissipative scale, but overd amped 
relat ive to the 'integral' (driving) scale, we can apply the Kolmogorov cascade principle (jFrisch 
19971 ). that motion in the inertial range is independent of the mechanism of dissipation. It is 
therefore determined by the rate of dissipation £ but it does not depend on ly. The moments of 
the relative velocity there fore depend only u pon £ and 7. For the second moment, dimensional 
considerations imply that (jMehlig et a/.l 120071 ) 

£ 



K- 



7 



(5) 



where Kolmogorov's 1941 theory of turbulence suggests that K a universal constant, b ut in practic e 
K should have a weak dependence upon Reynolds number due to intermittency effects (jFrischlll997l ). 
The constant K can be determined by simulation of Navier-Stokes turbulence, but this has not yet 
been done. We use if = 1 in the remainder of this article. 



Volk et al\ (|l98d ) (see also (jMizuno et a/.lll988l ). iMarkiewicz et ah (|l99ll )) have discussed rela- 



tive velocities in the large Stokes number limit, expressing their results in terms of properties of the 
integral velocity and timescales: their results are equivalent to dSJ, but are harder to apply because 
the precise form of the integral scale motion is highly uncertain, whereas the rate of dissipation in 
([5]) is determined directly from the accretion rate and the mass per unit area in the accretion disk. 



4. Properties of dust-particle aggregates 
4.1. Damping rate 

During the kinetic phase of planet formation, the size a of the suspended particles is usually 
smaller than the mean free path of the gas, A. In this case the drag force is proportional to 
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PgCsA{u — r), where pg is the gas density, A is the an gular average o f the projected area of the 
particle and Cg is the mean molecular speed of the gas (jEpstein Ill924l ). In this case the damping 
rate is 

PgCsA 



7 



m 



wher e m is the ma ss of the particle. In the case of spherical particles, the coefficient in eq. 
3/4 (jEpstein Ill923 ^. 



(6) 



IS 



4.2. Structure of aggregates 



The suspended particles are expected to consist of dust particles, of typical size oq, which are 
weakly bound together by electrostatic and van der Waals interactions. These aggregates may have 
a fractal structure, such that the number of particles in an aggregate of linear dimension a is 



N 



D 



(7) 



where -D is a fractal dimension. The damping rate for an aggregate of N dust grains with size ao 
composed of material of density pp is therefore estimated as 



7 



PgCs 
PpO-O 



N— 



70 



ao 



2-D 



(8) 



Since St = 1/(7t) we have St oc (a/ao)^ ^ and conclude that if D < 2, the Stokes number does 
not increase as the mass of the dust aggregate increases. 

In the early stages of aggregation, we expect that the mean free path of the suspended particles 
is very large compared to their size, and the relevant fractal growth mechanism is ballistic aggrega- 
tion (where particles incident with random linear trajectories adhere to the aggregate on contact) 
rather than diffusion-li mited aggregation (whe re an incoming particle executes a random walk). 
Numerical experiments (jNakamura et a/.lll994l ) on ballistic limited aggregation suggest D ~ 2.98 
for a 'ballistic particle cluster aggregate' (BPCA) model (particles adhering to a stationary aggre- 
gate) and D ~ 1.93 for a 'ballistic cluster cluster aggregate' (BCCA) model (where aggregates of 
particles move and collide with each other). (Note that the use of the term 'cluster' in the desig- 
nations of the BPCA and BCCA differs from that used in this paper, as explained at the end of 
section [T]) . 



However, experiments on agitated dust grains (jBlum et a/.l l200d : iKrause &: Blumll2004l ) have 
shown evidence of aggregates with a lower fractal dimension, D « 1.4. The lower fractal dimension 
observed in these experiments may be related to electrostatic effects: for example if agitation of 
dust particles causes them to acquire dipolar charges, they will tend to form strings (dimension 
D = 1), as illustrated in figure [2j The timescales on which these aggregates are observed is very 
short compared to the astrophysical context which we consider (and also the number of particles 
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in the aggregates are not very large). The appropriate value for D is, therefore, highly uncertain. 
It is probably not universal throughout the accretion disk structure. 

Our interest is in aggregates which undergo energetic collisions (which are sufficient to cause 
fragmentation in some cases), and which have a long time to relax into a compact configuration with 
lowest energy. In the following we therefore usually assume that the aggregates are predominantly 
in compact configurations, with Z) ~ 3, but in many cases below we continue to give formulae for 
general values of D. 



4.3. Strength of aggregates 

The other distinctive aspect of the dust-grain aggr egates concerns the i r frag il ity under colli 



sions. There is a substantial literature on this topic: see lDominik &: Tielena (|l997l ). lYoudinl (j2003l ) 
and references therein. In the following we introduce a new treatment, valid in the limit where the 
grains are very small. 

We mentioned earlier that the aggregates are bound together by van der Waals forces, which 
means that they must be very fragile. For a collision between two aggregates of size ai and 02, there 
will be a velocity Vqt above which the average change in the mass of the larger particle is negative. 
It is extremely hard to estimate this critical velocity for collisions of aggregates. We estimate this 
critical velocity below for the collision of two grains. We will assume that the critical velocity is 
approximately independent of size, so that the two-grain estimate is sufficient for the general case. 
In the following we therefore use Vcr to refer to the mean value of the relative velocity below which 
pairs of grains remain bound. 

We motivate this assumption by the following argument. When two aggregates collide, initially 
the collision only influences the grains on the surface of the aggregates in the vicinity of the point 
of collision. If the aggregates collide with a relative velocity v, some of the grains will acquire 
relative velocities which are comparable to v. If v/vcr is significantly greater than unity, grains on 
the periphery of the colliding aggregates will be able to escape. 

In order to estimate the critical relative collision velocity above which fragmentation occurs we 
therefore estimate the critical velocity Vcr above which two dust grains do not remain bound after 
a binary collision. This will be estimated by equating the kinetic energy of the collision process 




Fig. 2. — Dust grains in a low gravity environment are found to form aggregates with a low fractal 
dimension. This may be due to particles having an electric dipole, causing them to form chains. 
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between two dust grains, treated as spheres of radius oq, to the binding energy. 



The binding energy is us ually es t imate d by a theory proposed by IChokshi et al\ (jl993l ) (and 
discussed in simpler terms by lYoudinl (|2003l )). which assumes that the spheres deform so that they 
are in contact on a small circular patch. By balancing the force required to deform the spheres 
against a derivative of the surface energy released by bringing them into contact, the equilibrium 
deformation and hence the binding energy are estimated in terms of a Young's modulus for the 
spheres, Y , and their surface energy of attraction, a. Omitting the dimensionless prefactors (which 
depend upon making more precise definitions of y , a) the surface deformation 6 and binding energy 
£'b are estimated to be 



y2 



1/3 



5„4\ 1/3 



a a. 



y2 



(9) 



The latter expression gives an estimate for the critical velocity Vcr, obtained in lChokshi et al\ (|l993l ): 



1/6 



(10) 



Using val ues for quartz of Y = 5.4 x 10^°Nm~^, cr = 2.5 x 10^^ Nm^^, pp = 2.6 x lO^kgm"^ 
quoted in IChokshi et al\ ( 19931 ). we estimate 5 ~ 3 x lO^^^m and Vcv ~ 0.2ms"^ for grains of 
size Co = 10~^m. The validity of this approach is questionable when the value of the deformation 
is small compared to the size of an atom, as is found in this example. We therefore propose an 
alternative estimate of the binding energy. 

We assume that there is a van der Waals binding energy 5Ei between atoms, which operates 
over a range /, which we equate with the typical interatomic distance. We assume that there is a 
binding energy SEi/t^ per unit area for all parts of two spheres which are within a distance 21 <^ 
of each other, where 6E is the van der Waals binding energy between two atoms, and I is the range 
over which it acts. Two spheres in contact have a surface area nlao with separation less than 21 (see 
figure [3|), so that the binding energy is £'b ~ nSEao/l. Unless the grains are made of very plastic 
material, particles are not expected to aggregate if their (centre of mass frame) kinetic energy is 
large compared to the binding energy. Writing E^ ~ mw^ ~ PpCLQi'cn 'we estimate a critical velocity 



5Ei/p^la, 



(11) 



Taking 5Ei 10^^^ J, Z = 3 x 10^^'^m and ao = 10^'m, we find that ~ 3 x lO^^ms^^ for binding 
of microscopic grains. Thus collisions between aggregates with relative velocities of the order of 
30cms~^ have the potential to dislodge particles from an aggregate. 

It is instructive to compare the microscopic and the macroscopic models, equations (jlOp and 
TTI ). The surface energy and Young's modulus can be related to microscopic parameters by writing 



a ~ SEi/P, where 6Ei ~ 10 J is the van der Waals interaction energy and Y ~ 6E2/l^, where 
6E2 ~ 10"^^ J is the characteristic energy scale for covalent bonds: these expressions give values 



Fig. 3. — Two spheres in contact have a surface area vr/ao with separation less than 21. 



which are consistent with the values for a and Y for quartz quoted above. Expressing equation 
pUj) in terms of microscopic quantities gives 

1/3 



^E.Mr^] (f) ^ (12) 



We see that the two estimates (|10p and (jlip differ by dimensionless ratios with small exponents, 
and the values obtained from the two formulae are therefore typically quite similar. Whichever of 
equation (fTTI) or (flO]) gives the larger value should be used to estimate the critical velocity. The 
material parameters can make a significant difference, in particular ice has a much larger surface 
energy and is softer tha n quartz, and the pr edicted critical velocity is approximately an order of 



magnitude larger for ice (jChokshi et a/.lll993l ). 



There are significant sources of uncertainty in these estimates for the critical velocity for 
binding single grains, which are very hard to quantify. One is the effect of surface roughness. This 
will reduce the effective contact area between two grains, and therefore reduce Vcr- The effect of 
electrostatic charging of dust grains and aggregates is harder to quantify. It is possible that the 
microscopic grains will all tend to acquire a positive charge, due to ionisation by energetic photons. 
This would reduce the binding energy. Another possibility is that grains acquire random charges 
through frictional charge transfer in collisions. This may increase the critical velocity for pairs 
of grains, because of electrostatic attraction of oppositely charged particles. This effect will be 
reduced by cancellation of charges for larger aggregates. 

Experiments to determine the critical velocity fo r collisions of grains have been performed 



by projecting small silica spheres at a glass surface (iPoppe &: BlumI 119971 ). These experiments 
have usually showed that the spheres adhere to the surface at velocities which are approximately 
ten times higher than those predicted by equation (jlOp . The reasons for this discrepancy are not 
understood. Although the effects of charges on the particles were quantified in these experiments, 
it is possible that effects of electrostatic charges on the surface might account for the discrepancy. 



5. Turbulence in protoplanetary disks 



We now consider estimates relating to dust particles suspended in gas surrounding a growing 
star. Observational evidence suggests that there is usually sufficient angular momentum that the 



- 12 - 



gas fo rms an accretion disk and we use the steady-state disk model described bv lShakura &: Sunvaev 



(jl973l ). The kinetic and potential energy of the material in this cloud must be dissipated into heat 
and radiated away in order to allow the material to fall into the growing star: we consider the 
case where the disk is primarily heated by dissipation, rather than by radiation from the star. It is 
known that stars can form quite rapidly (over a timescale of tacc ~ lO^yr). The process is so fast 
that either shocks or turbulent processes must play a role in the dissipation of energy: the rate 
of dissipation by laminar flow is too small by many orders of magnitude. It has been argued tha t 



magneto-hydrodynamic mechanisms provide the necessary large-scale instability (lArmitagd 120071 ) . 
but the gas is only weakly ionised and the small scale dissipation is expected to be described by 
conventional hydrodynamics, for which the Kolmogorov theory of turbulence is applicable. Shocks 
dissipate a finite fraction of the kinetic energy of a gas almost instantaneously, so that their role 
must be short lived and throughout most of the accretion process turbulence is the dissipation 
mechanism. 

In the following we consider a star of mass M (taken to be one solar mass) surrounded by a 
protoplanetary nebula, which is collapsing into the star at a rate M. We assume 

M = lO^'^AMyr^^ . (13) 

Most discussions of protoplanetary accretion systems assume that A ~ 1. We assume that the 
accretion zone is predominantly molecular hydrogen (we take the mean molecular mass of the gas 
molecules = //mn, with n = 2.34 and mn = 1.67 x 10~^^kg being the mass of atomic hydrogen). 
For estimating the mean-free path of the gas, we assume that the molecular collision cross section 
is 5coi = 2 X lO^^^m^. These and other assumed parameter values are collected in tabled! We use 

cl = kBT/rrig (14) 

to relate Cg, a characteristic speed which is of the order of the speed of sound, to temperature. 

The accretion disk is characterised by its mass density per unit area 5](i?) at a distance R 
from the star, its scale height H{R), its rate of dissipation per unit area Q{R), and its surface 

Table 1. Assumed values for parameters. 



Variable Symbol Value 



Mass of gas molecules 


rUg = fimn 


fi = 2.34 


Molecular collision cross section 


Scol 


2 X lO^^^m^ 


Mass of the star 


M 


1.99 X 103°k^ 


Typical size of interstellar dust grains 


Go 


lO-'^m 


Density of suspended particles 


Pp 


2 X lO^kg/m' 


Mass fraction of suspended particles to gas molecules 


K 


3/170 
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temperature T{R). Using the continuity equations for flow of mass and angular momentum within 
a thin accretion disk which is in a quasi-stationary state, the r ate of dissipation per unit area Q at 



radius R can be obtained in terms of the rate of accretion, M (jShakura &: Sunvaevlll973l ): 

3 GMM 



Q{R) 



' -R 



(15) 



Here Rc is the radius of the core of the accretion system. In the following we concentrate on the 
region R/ Rc ^ 1, and do not carry forward the final factor of the above expression. Note that 
this formula is independent of the actual mechanism of dissipation. The rate of dissipation per 
unit mass is obtained from the mass per unit area S(i?) of the accretion disk at radius R, and the 
density of gas Pg{R) at the mid-plane of the disk is related to the scale height: 

where in the second expression we assume a Gaussian density profile with variance H. Con- 
sider next how to estimate the height of the disk, H{R). The gas is in Maxwell-Boltzmann 
equilibrium in the gravitational potential and we assume that the contribution to the potential 
from the disk itself is negligible. The potential energy at a distance z from the mid plane is 
<I>(z, i?) = GMnigz'^ /2R^ , so that the gas density at distance z from the mid plane of the disk is 
p{z, R) = Pg{R) exTp[—z'^/2H^{R)], the scale height of the disk being 



where i^{R) = \J GMj R^ is the Keplerian orbital angular frequency for a mass M at radius R. 

We assume that the surface temperature T{R) is determined by the rate at which the ther- 
mal energy created by dissipation of turbulent motion can be radiated away. This assumption is 
justifiable as long as A in equation (fT3]) is not too s mall, but when A <C 1 the rate of heating by 



radiation from the star is expected to be significant (| Armitagd 120071 ) . Assuming that the surface 



of the disk behaves as a black body with emissivity e, one obtains 

Q{R) = 2eaT\R) ~ (18) 

where a is the Stefan-Boltzmann constant. Assuming e ~ 1, one obtains the radial temperature 
profile 

where Rq is a convenient reference radius, which we take to be Rq = lAU = 1.5 x 10^^ m, and 
T{Rq) is the temperature at Rq with A set equal to unity. Given the radial dependence T[R) of 
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the temperature, the radial dependence of the velocity of sound and of the disk thickness can be 
determined. Using equations (jlSp and (|17p we find 



.(i?) = c.(i?o)Ai/«(|^) 



H{R) = H{Ro)A'/'(-^] 



-3/8 



9/8 



(20) 



(21) 



(Again, the dependence on A is not shown explicitly in the list of arguments.) In order to estimate 
the radial dependence of other quantities it is necessary to determine the radial dependence of the 
density. This is achieved as follows. Applying the continuity equation for angular momentum leads 
to a relation between E and an effective kinematic viscosity (termed the eddy viscosity), fefr, that 
is th e diffusion coefficient chara cterising the transport of angular momentum across the accretion 
disk (jShakura &: Sunvaevlll973l ): 

M = 3Tri^cs{R)T.{R) . (22) 



The eddy viscosity in turn is determined by the size L of the largest eddies in the turbulent flow 
(the integral scale). We estimate ~ L'^/ti, where t_L is the eddy turnover time for an eddy of 
size L. The Kolmogorov scaling argument implies that ~ (L^/iS)^/'^, so that we estimate 



It is natural to assume that L is of the order of H{R) and we write 

L = £H{R) 



where £ is a parameter. Taking (|T7j) . (|T8|) . (p2|) and (j23|) together, we find 



2V2 M 



(23) 



(24) 



(25) 



97r£2 

Usually it is assumed (jShakura &; Sunvaevlll973l ) that u^s ~ ac^H (where a < 1 is a dimensionless 



coefficient). Our own approach is equivalent to this 'a-prescription': using ([251) in equation ([22 
we obtain 



2^/2 



(26) 



Thus we see that the a-presription is equivalent to assuming that the integral scale of the turbulence 
is smaller than H hy & fact or i k, ^fg. It is wide ly accepted that may many observations are 
compatible with a ~ 10^^ (IHartmann et a/.lll998l ). corresponding to L ~ i^/lO. Accordingly, 
wherever we quote numerical values for quantities without specifying how they scale with £, we 
have set I = 0.1. An advantage of our '^-prescription' is that it makes the physical nature of 
the adjustable parameter clearer than for the standard 'a-prescription'. The rather small values 
of a indicated by observations could be indicative of a fundamental problem with a theory. Our 
alternative approach is reassuring because it indicates that a more physically transparent parameter, 
I = L/H, is not in fact a very small number. 
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These results allow us to determine the power-law radial dependence of other quantities. The 
results for those determining the gaseous component of the accretion disk are summarised in tabled 
writing a generic variable in the form 

X = X{Ro) ( — ] A^'^^ . (27) 

where Rq = 1 AU = 1.5 x 10^^ m, and X(Rq) is the value of X at radius Rq with A = i = 1. 



6. Relative velocities, collision rates and Stokes trapping 

Next we consider the behaviour of dust grains suspended in the gas forming the accretion 
disk. We consider the case where the dust grains are almost all sub-micron sized particles: we 
assume they are composed of material with density Pp = 2 x lOfkgm"^, and that these are initially 



spherical particles of radius ao = 10 ^m. Following iHavashil ( 198ll ). we assume that the mass ratio 
K of suspended particles to gas molecules is k = 3/170 ^ 0.018. 

The results described above can be used to estimate collision rates. We consider two cases. 
First, we estimate the collision rate when the Stokes number is small , and when the collision 



mechanism is Brownian diffusion (we shall see that advective collisions ([Saffman Turneij 119561 ) 
make a negligible contribution). Second, we also require the rate of collision for large Stokes number. 
Both estimates require the number density of dust grains. We write 

n^^(^y"=no(^y" (28) 
mo \aoJ \aoJ 

where D is a fractal dimension and k « 0.018 is the ratio of the mass densities of condensed to 
gaseous matter and tuq = Airppa^/S is the mass of a microscopic dust grains; the second equality 
defines no- 

The mean speed of particles due to Brownian motion is 

where the second equality defines vo- Our estimate of the collision rate due to Brownian diffusion 
for a general value of the fractal dimension of the aggregates is 

4-3D 



TZd ~ 47r\/2na^(fd) ~ AirV^noalvo ( — ) 



2 



4-3D 

a ^ 



and the factor TZ^o is given in table [3l The collision rate due to advective shearing is approximately 

7^a . (31) 

r 



- 16 - 



The gas-kinetic collision rate (taking K = 1 in equation ([5])) is 





2-D 
2 




) 



2-D 
2 



(32) 



The Stokes number is 



1 




) 



D 



2 



St 



(33) 



7r 



All of these quantities have a power-law dependence upon radius and particle size: for the generic 
quantity X we write 



(again, it is to be understood that X{Rq, oq) is evaluated for A = £ = 1). The quantities determining 
the collision rates are collected in table [3j 

The growth of the aggregates due to collisions could be be limited by two factors. Firstly, the 
rate of collisions usually decreases as the size of the aggregates increases and the collision rate could 
become negligible when the particles grow to a sufficiently large size. Secondly, large aggregates 
could be more vulnerable to being fragmented upon collision. 

To illuminate a discussion of whether these limitations occur in practice, we now consider the 
properties of particles with size a chosen so that St = 1. For compact particles {D = 3) we find 
that their size is 



with a*(i2o) ~ 6 X 10~^m: for I = 0.1 this leads to particles of size a* ^ 25 /um at i? = 10 AU. It 
is also of interest to evaluate the collision rate of these particles. Below we describe a mechanism 
which will lead to the aggregates becoming compact, so that we concentrate on the predictions of 
this calculation when D = 2>. The approximate values are listed in table [H Note that the diffusive 
collision rate is very small, but it increases abruptly (by approximately five orders of magnitude, 
for the case where £ = 0.1) when the size of the particles increases such that the Stokes number 
exceeds unity (see figure H)^ a)). Also, note that the relative velocity for collisions of particles with 
Stokes number unity is comparable to the critical velocity for microscopic dust aggregates to be 
fragmented. 

For compact aggregates (D w 3), the collision rate for the gas-kinetic mechanism is the much 
larger when St ^ 1, and it also decreases less rapidly as the size of the particles increases. If the 
lifetime of the accretion system is tacc ~ lO^yr, the growth process may be assumed to be ended 
when T^gtacc = 1- The estimate for TZ^ above indicates that the gas- kinetic collision mechanism is 
sufficiently rapid that it does not limit the growth of aggregates, at least until they are so large that 




(34) 




3/8 



(35) 



-17- 



the approximations in equation (|33p . such as assuming that the drag is determined by the Epstein 
formula, break down. We conclude that kinetic factors do not limit the growth of aggregates, once 
particles become sufficiently large that St > 1. Note, however that for compact particles at the 
size where St = 1, the diffusive collision rate is extremely slow (of the order of one collision per 
10^ years when £ = 0.1). The slowing down of the collision rate as the particle size increases is 
therefore quite close to becoming an insurmountable bottleneck. 

As the size a of the aggregates increases, so does their relative velocity (Au), for example 
compact aggregates of size a = 10 cm collide with a relative velocity of approximately 10ms~^ if 
we take i = 0.1. Such a collision would certainly cause the aggregates to fragment upon impact. 

We have arrived at the following picture of the initial stages of the growth of dust particle 
aggregates. Initially, dust particles will aggregate by Brownian diffusion. This process continues 
until the Stokes number of the larger particles is of order unity. Then the velocity of these particles 
starts to separate from that of the gas due to the formation of fold caustics. This greatly increases 
the relative velocity between these particles and smaller ones in their vicinity (figure Hl^a)): the 
larger particles will rapidly 'sweep up' the remaining small particles in their vicinity. However, 
when these larger particles collide with each other they have sufficient kinetic energy that they 
cause each other to fragment, so that particles with large Stokes number fragment to replenish the 
supply of particles with small Stokes number (figure|H|b)). This results in a steady-state distribution 
of particle sizes, corresponding to Stokes numbers of order unity (illustrated in figure [D^c)): we call 
this the 'Stokes trap'. 

7. Mechanism for avoiding the Stokes trap 

In the following we propose four mechanisms by which the Stokes trap might be avoided. 



7.1. Escaping the Stokes trap by clustering 



It is known that small particles suspended in a turbulent flow can cluster, even if the fl uid flow 
i s inco mpressible, due to effects of the inertia of the particles. This effect was proposed by iMaxey 
(jl987l ). who suggested that heavy particles would be 'centrifuged' out of vortices and would tend 
to cluster in regions of low vorticity: this effect is often referred to as 'preferential concentration' 
because it is thought that the particles cluster in regions of low vorticity. 

The argument by Maxey does not give a quantitative description of the nature of the clustering 
effect at long times. ISommerer &: OttI (jl993l ) suggested that particles in random fluid flows cluster 
onto a fractal se t of di mension Dd and proposed using the Lyapunov dimension (first defined by 
Kaplan &: Yorkd (119791 )) to characterise this set. Recent numerical experiments (iBec et a/.ll200a ) 
have studied how the (Lyapunov) fractal dimension Dl of particles in a three-dimensional turbulent 
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10"'' lO"'" 10"^ 10"' 



Fig. 4. — Schematic illustration of the Stokes trap, (a): Relative collision velocity of equal-sized 
aggregates versus their size a. (b): Small aggregates grow by collisions, but larger aggregates tend 
to fragment upon collision, (c): Corresponding stable particle-size distribution. 



flow varies as a function of Stokes number: the minimum value is -Dj. ~ 2.6 for St = .55, -Dl = 3 
for St > Stc ~ 1.7 and Dl — > 3 as St ^ 0. Theoretical work (jPuncan et a/.ll2005l ) has shown 



that the fractal clustering may be explained without the centrifuge (or 'preferential concentration') 
effect: a model for which it is absent gives a minimum dimension Dl ~ 2.65, in good agreement 
with simulations. 

It has previou sly been suggested that this clustering effect may play a role in planet formation 



(|Cuzzi et a/.ll200ll ). but there is a difficulty. The clustering effect only occurs when St = 0(1), 
but most aggregation processes produce a broad (and time-dependent) distribution of sizes. The 
clustering mechanism would then only act on a small proportion of the particles. Our Stokes- 
trapping mechanism avoids this difficulty. 

The Stokes-trapped dust aggregates will be clustered, but if the clustering is confined to length 
scales below the Kolmogorov length r], the effect is not sufficiently strong to initiate gravitational 
collapse: the mass of solid matter within a Kolmogorov length is 

51/16 



M 



1.0 X IQ-^kgA-^^/^^i 



E_ 

Ro 



(36) 



The size of the clusters created by the 'preferential concentration' effect cannot significantly exceed 
this value, and this mass is much too small for gravitational effects to become significant. 



We thus conclude that although the clustering effect is promoted by the Stokes trap, clustering 
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of dust aggregates onto a fractal set cannot provide inhomogeneities which are sufficient to trigger 
gravitational instability in protoplanetary systems. However, it is possible that this clustering effect 
could influence the scattering of electromagnetic radiation by the dust grains. 



So me other mechanisms for clustering have been proposed, discussed in iBarge &: Sommeria 



(jl995l ). Ijohansen et a\ (120071 ). We argue that these mechanisms can be discounted because they 
assume the existence of much larger aggregates than the Stokes-trapped size. 



7.2. Escaping the Stokes trap by chemical processes 

We propose that the mid-plane of the accretion disk can be at a much higher temperature 
than its surface. The dust aggregates circulate between regions with different temperatures. In 
the lower temperature regions, water and organic compounds can condense onto the aggregates. 
When the particles circulate to higher temperature regions chemical reactions may occur. This is 
expected to both compactify the aggregate and to make it much more resistant to being fragmented 
in collisions. It is also possible that the temperature in regions close to the star is sufficiently high 
that at least some of their component materials are melted. This will also produce much more 
robust and compact aggregates, able to withstand collisions at much higher relative velocities. 

We have assume that the temperature of the accretion disk is determined by dissipative heat- 
ing, and we have estimated the surface temperatures, which are quite low. However, if the disk is 
optically thick in the spectral region corresponding to the Wien wavelength, the mid-plane tem- 
perature may be considerably higher than the surface temperature. 

The optical thickness W is the ratio of the height of the disk to the photon mean free path. 
The optical properties are dominated by the dust particles \IW ^ 1, and when the particles are 
larger than the wavelength we may assume that their optical cross-section Sopt is approximately 
the square of their linear dimension: 5opt ~ a^. If the density of particles is n, the optical mean- free 
path is Aopt ~ l/(""S'opt)- Thus for compact particles {D = 3), the optical thickness is 

W ~ Hna^ . (37) 

The ratio of the interior temperature Tint of the accretion disk to its surface temperature T is 

^ ~ W^'^ . (38) 

Very small particles, with dimension small compared to the Wien wavelength corresponding to 
temperature T, do not absorb or scatter radiation effectively. If the aggregates are very large, 
their number is correspondingly reduced, so that they make a smaller contribution to the optical 
thickness. From the data above, we see that the Stokes-trapped particles have sizes which are 
comparable to, but somewhat larger than, the Wien wavelength corresponding to the surface tem- 
perature. The Stokes-trapped particles are therefore of approximately optimal size to increase the 
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optical thickness. The optical thickness corresponding to Stokes-trapped particles is therefore (for 
compact particles, D = 3): 

W* ~ ~ 410 ( ^) . (39) 

a* \Ro J 

The temperature at the centre of the disk is therefore expected to exceed the surface temperature 
by a factor of approximately four in the range R = 1 — lOAU. 

The chondrules, grains of typically submillimetre size which are found in many meteorites, 
show evidence of having been heated to very high temperatures, causing melting. The temperatures 
required are rather higher than those which we predict above (except very close to the star). Also, 
the exposure of chondrules to high temperatures may have been brief compared to the circulation 
time tL, otherwise the chondrules may have evaporated. 



7.3. Escaping the Stokes trap by forming tenuous aggregates 

We have already remarked that if the fractal dimension D of the dust aggregates is D < 2, 
then the Stokes number does not increase as the size of the aggregate increases. Such aggregates 
would avoid the 'Stokes trap', but there are three reasons why these tenuous aggregates might not 
provide a satisfactory solution to the problem of planet formation. 

First, it might be thought that as such aggregates grow their collision rate would become so 
small that growth would, for all practical purposes, cease. In section |3] it was determined how the 
collision rate for fractal aggregates with dimension D scales as a function of their characteristic 
size, a: we found TZd cx a*^^~'^^^/^. The collision rate therefore increases with the aggregate size 
when D < 4/3, because such very tenuous aggregates are so extended that they entangle each 
other. Experiments on growing aggregate s in negligible grayity have produced ob jects with very 



low fractal dimensions of roughly D = lA (iBlum et a/.ll2000l : iKrause &: Blumll2004l ). rather smaller 



than the values predicted for kinetic aggregation processes. The dust aggregates observed in low- 
gravity experiments are only observed for a short period compared to the time between collisions 
between dust aggregates, so that the structure of the aggregates in protostellar accretion systems 
is uncertain, but there is probably no kinetic restriction on the growth of very tenuous clusters. 

A second possibility is that the tenuous aggregates are torn apart by the effect of the shearing 
motion of the fluid. We estimate this effect as follows. Consider the force required to split an 
aggregate of fractal dimension D, size a, composed of particles of size oq. The number of bonds at 
an equatorial plane of the aggregate is approximately 

Each bond has a breaking strain Fb. Grains in contact have an energy per unit area 5Ei/P per 
unit area, acting over a range I. The binding energy is E'b = TT5EiaQ/l (see section . Writing 
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Eh = -^b^ we estimate 

Fh ~ . (41) 

The drag force on a single grain due to fluid motion with a speed u relative to the fluid is F = rriQ^ou. 
When D < 2, the drag force pulling the aggregate apart is comparable to the sum of the forces 
which are predicted by applying this formula to each grain. The relative velocity of the gas at 
opposite sides of the aggregate is a/r, so that the force acting to break up the aggregate is 



T \aoJ 

The ratio of the applied force to the critical force is Fbr/iVb^b! note that this ratio is independent 
of D (provided D < 2). Setting this ratio equal to unity we estimate the size Omax of aggregates 
which will be torn apart by shearing forces: 

a„,ax«78mA-9/i6W— j . (43) 

The estimates leading to this result are valid provided the aggregate is small compared to the 
Kolmogorov length: a <C r/: this condition is only marginally satisfied. It is not clear how the 
shear forces increase with aggregate size once the aggregate is large compared to the Kolmogorov 
correlation scale, but it is hard to escape the conclusion that tenuous aggregates are torn apart by 
shear forces when they become sufficiently large. The mass contained in these tenuous aggregates 
is clearly comparable to that estimated in equation (j36p and is too small to initiate gravitational 
collapse directly. 

A third problem is that aggregates with Stokes number less than unity are always advected 
with the gas. Even if they could become very massive, they would never collapse to the mid-plane 
of the disk, as long as the motion of the gas remains turbulent. 



7.4. Escaping the Stokes trap by reducing the accretion rate 

If the accretion rate M is reduced (by choosing a small value for the parameter A in equation 
(jlSP ). there is a corresponding reduction in the turbulence intensity £. We might expect that there 
would be also be a reduction in the relative velocities of suspended particles. This latter expectation 
is false, because a reduction in the rate of accretion results in a reduced density of the gas, and 
the suspended particles are more lightly damped by the gas (note that if 7 decreases the relative 
velocity increases: see equation ([5])). 

This principle is illustrated by results in table HI which show that the Stokes-trapped particle 
size decreases as the accretion rate parameter A decreases (because 5/\_ > for a*). A further 
demonstration comes from considering the size of particles which are not fragmented on collision. 
If the relative velocity is less than the critical velocity Vcr for disintegration of dust aggregates, then 
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they continue to grow on collision. The size of particles for which the turbulent relative velocity 
exceeds Vct is 

Ocr = 1.0 X 10"^ s^m-^vl _ (44) 

We see that, contrary to intuition, the attainable particle size decreases as the rate of accretion 
decreases. 



8. Settling and gravitational collapse 

Even if the Stokes trap is circumvented, the particles must grow to a substantial size in order 
for them to settle to the mid-plane of the accretion disk. In the following we illustrate this fact 
by estimating the size at which particles start to settle to the mid-plane, according to our model. 
We find that settling only starts at large particle sizes (approximately 10 cm, where the use of the 
Epstein formula for the damping rate becomes questionable). We also estimate the critical height 
of the dust layer for the onset of gravitational collapse, and find that this is a small fraction of the 
height of the gas layer. 



8.1. Critical size for settling 



The integral scale fluctuations with velocity ul will advect aggregates away from the mid plane. 
This effect is opposed by gravitational attraction to the mid plane. The gravitational acceleration 
at a distance z from the mid plane is 



g{R,z) 



GM 



(45) 



If the integral timescale of the turbulence, is large compared the the relaxation time 7~^, then 
the aggregates only settle to the mid-plane if their terminal speed ut = g{R, z)/^ exceeds ul- There 
is thus no settling at all unless 7 < 7s, where 



7s 



UL 



1.8x10-.-'.- (I) 



-3/2 



(46) 



In the case of compact aggregates {D 
plane when their size exceeds 



3), we conclude that particles start to settle to the mid- 



6.0 X 10~^m 



R 

Rq 



-3/4 
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When considering very large aggregates, the damping rate 7 may be sufficiently small that "yti <C 1- 
In this cas e, the effect of the turbulence ra ust be modelled as a sequence of random impulses on the 
aggregate (jCuzzi et gj 119931 : lYoudinl 120071 ). We argue that, because of the Stokes trapping effect, 
the aggregates never grow to a sufficient size for this alternative approach to become relevant. 
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8.2. Critical height for gravitational collapse 

A stationary gas of solid particles with mass density psoi is unstable against gravitational 
collapse on a timescale tcoU ~ (G/Ogoi)"^^^- If the gas consists of particles with typical relative speed 
Av, then gravitational collapse (the Jeans instability) occurs on length scales greater than Avtcon- 
In the case of an accretion disk, the question of gravitational instability can be complicated by the 
rotational motion and the finite thickness of the disk. 

The gravitational stability of a uniform thin disk of area density T, may be described by giving 
a dispersion relation for the frequency w of a density perturbation with wavenumber k 

co^ = Q'^- 27rGS|A:| + Vsk"^ (48) 



where Vs is a two-dimensional sound velocity (jBinney &; Tremaind Il988l ) . The system becomes 
gravitationally unstable when u>'^ < for any value of k. For a collisionless system, we have 
Vs = 0, and the instability arises when the term containing the gravitational constant is larger in 
magnitude that the centrifugal term. The instability is favoured by choosing a large value for k, 
but if we consider very large values of k, the approximation of treating the mass distribution as 
two-dimensional fails. We therefore assume that the largest possible value of k is 2it/H, where H 
is the height of the disk. 

The solid material in the disk (with area density kT,) therefore becomes gravitationally unstable 
when its scale height, Hsoi is less than a critical value which is given (approximately) by 

= - 3.3 X lO^m ( ^) . (49) 



■ \Rq) 

Note that Her/ H ^ 1; implying that to trigger a gravitational instability, particles must be very 
much heavier than particles of the size Ocr (given by equation (j47|) ). 



9. Concluding remarks 

In this paper we have analysed the consequences of turbulence for a standard model of the 
formation of planets by aggregation of sub-micron sized dust. This led us to introduce the concept 
of the Stokes trap, a mechanism which limits the growth of dust aggregates in turbulent protostellar 
accretion disks and thus constitutes a barrier to planet formation. The particle sizes predicted by 
the Stokes trap mechanism are so small that turbulent fluctuations would never allow the particles 
to settle to the mid-plane of the disk and achieve sufficient density to trigger gravitational instability. 
The Stokes trap must therefore be avoided if the planets are to form. 

We have proposed four possible ways of avoiding the Stokes trap and assessed these mechanisms 
in the light of Kolomogorov's scaling principle for turbulence in combination with a standard 
model for the accretion disk. Particles clustering due to 'preferential concentration' cannot produce 
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sufficiently large clusters. Chemical processes which consolidate aggregates will strengthen them, 
but not to the extent that large aggregates would be invulnerable to fragmentation by collision. If 
the aggregates arc very tenuous fractals, with dimension D < 2, their Stokes number remains very 
small, but the aggregates may be torn apart by shearing forces when they become sufficiently large, 
and they will not settle to the mid-plane of the disk. Finally, we observed that as the disk becomes 
less active (that is, as M decreases) the relative velocity of the suspended aggregates increases, so 
that the difficulty is not removed as the disk becomes quiescent. 

Even if the aggregates escape from the Stokes trap and continue growing, our estimates indicate 
that they must reach very large sizes before they can overcome the effects of turbulence and start 
to settle to the mid-plane. Our estimates indicate that the turbulence intensity would have to be 
reduced by many orders of magnitude before this could happen. 

We conclude that a model for planet growing by aggregation of small dust grains is highly 
problematic. It is therefore necessary to examine other possibilities. The difficulties with the theory 
are avoided if the protoplanetary accretion disk already contains large objects with a gravitational 
escape velocity which is large enough to prevent them from being disrupted by collisions with dust 
grains and aggregates. A star forms when a cloud of interstellar gas satisfies conditions which allow 
it to undergo gravitational collapse. It is likely that at the same time as the star forms, other 
regions of the cloud in the vicinity of the nascent star will condense by gravitational attraction. 
Condensed objects much smaller than the star itself may be drawn into the accretion disk of the 
star, and eventually become planets lying on roughly circular orbits in the equatorial plane of the 
star. Thus we propose an alternative picture for the formation of planets which we term 'concurrent 
collapse', in which the planets are formed as gravitationally bound systems at the same time as the 
star. During the life of the accretion disk they become coupled by friction and mass transfer to the 
circular motion of the accretion disk. While in the accretion disk their structure may be radically 
transformed by processes such as further accretion of material. Their mass might also be reduced 
by ablation or evaporation of lighter elements. 

One possible criticism of this scenario is that the density p of the gas cloud at temperature T 
which collapses to form the star will be associated with a mass scale, the Jeans mass. 



It may be argued that the Jeans mass is a lower limit to the mass of objects that could form by 
gravitational collapse. However, as gravitational collapse proceeds, the density of some regions will 
increase, thus lowering the Jeans mass estimate which applies for further collapse. Images of gas 
clouds suggest that they are typically non-uniform in their density, and that they have density 
fluctuations spanning a wide range of spatial scales. The distribution of sizes of dense objects 
produced by gravitational collapse from such a non-uniform initial condition might be expected to 
be very broad. 

For the time being our 'concurrent collapse' hypothesis must be supported by the implausibility 




(50) 
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of the standard dust aggregation model. Because the additional hypothesis only concerns conditions 
at early stages of the life of the stellar accretion system, its implications for the eventual structure 
of a solar system are unclear. 

Support from Vetenskapsradet and the platform 'Nanoparticles in an interactive environment' 
at Goteborg university are gratefully acknowledged. 
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Table 2. Estimates of quantities characterising the gaseous part of the accretion disk, quoted 

to two significant figures. 



Variable 


Symbol 


Equation 


X{Ro) 








Surface temperature 


T 


eq.CSI) 


130 K 


-3/4 


1/4 





Speed of sound 


Cs 


eq.dSni) 


670 ms~^ 


-3/8 


1/8 





Disk height 


H 


eq.dlll) 


3.4 X 10^ m 


9/8 


1/8 





Surface density 


S 




280 kg m"2 


-3/4 


3/4 


-2 


Gas density 




eq-dlH) 


3.3 X lO^^kgm"^ 


-15/8 


5/8 


-2 


Dissipation rate 


£ 


eq-dlH) 


0.11 m^s"^ 


-9/4 


1/4 


2 


Gas mean-free path 


A 




0.42 m 


15/8 


-5/8 


2 


Kinematic viscosity 


V 


V = Acs 


280 m^s"^ 


3/2 


-1/2 


2 


Kolmogorov length 


7] 


eq.® 


120 m 


27/16 


-7/16 


1 


Kolmogorov time 


T 


eq.© 


52 s 


15/8 


-3/8 





Kolmogorov velocity 


UK 




2.3ms-^ 


-3/16 


-1/16 


1 


Integral velocity 


UL 


UL = L)l/3 


710ms"^ 


-3/8 


1/8 


1 


Integral timescale 


tL 


tL = L/uL 


4.6 X 10^ s 


3/2 









Table 3. Estimates of quantities characterising the suspended particles (quoted to two significant figures). 



Variable 


Symbol 


Equation 


X{Rq, ao) 


Sr 








Sa 


Damping rate 


7 


eq.® 


O.lls-i 


-9/4 


3/4 


-2 


2 


- D 


Number density 


n 


eq.dlll) 


7.0 X lO'^m^^ 


-15/8 


5/8 


-2 




-D 


Difi'usive velocity 




eq.^ 


1.4 X lO'^ms"^ 


-3/8 


1/8 







D/2 


Difi'usive collision rate 




eq.dHO]) 


1.8 X 10""^ s"^ 


-9/4 


3/4 


-2 


(4- 


3D)/2 


Advective collision rate 




eq-dUD 


1.4 X 10~^^s~^ 


-15/4 


1 


-2 


2 


- D 


Relative velocity 




eq.®, K=l 


0.98 ms"i 





-1/4 


2 


(D 


-2)12 


Gas-kinetic collision rate 




eq.(l32D 


3.6 X 10-6 s-i 


-27/16 


3/16 


1 


(2- 


-D)/2 


Stokes number 


St 


eq.dMl) 


0.17 


3/8 


-3/8 


2 


D-2 



I 

CO 

o 

I 



4. Estimates for the compact particles {D = 3) with St = 1. The results are of the form (j27l) . 



Variable 


Symbol 


Equation 


X(Ro) 






Si 


Stokes-trapped particle size 


a 


eq.(l35]) 


5.7 X 10"'^ m 


-3/8 


3/8 


-2 


Diffusive collision rate in Stokes-trap 




eq.dSO]) 


2.3 X 10-^s-i 


-21/16 


3/16 


3 


Advective collision rate in Stokes trap 




eq-dSI]) 


2.5 X lO^^^s-i 


-39/8 


17/8 


-8 


Gas-kinetic collision rate in Stokes trap 




eq.(l32]) 


3.6 X 10^6 s-i 


-27/16 


3/16 


1 


Turbulent relative velocity at St = 1 


(At;)* 




2.3ms"i 


-3/16 


-1/16 


1 



